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Using a large flavor number expansion and a gap equation for the pion mass the 
chiral quark-meson model is solved at the lowest order in the fermion contributions. 
In the chiral limit the tricritical point (TCP) is determined analytically. The 
softening of the sigma particle is verified at this point and is further investigated 
for a physical pion mass in the neighbourhood of the critical endpoint (CEP). 



1. Introduction 

Recently the location of the CEP in the fi — T phase-diagram of QCD was 
obtained for Nf = 2 + 1 dynamical staggered quarks with physical masses 1 . 
One might investigate this problem using an effective low energy model also, 
with some nonperturbative technique. Its low energy nature might prevent 
us from giving the accurate location of the endpoint and the shape of the 
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phase boundary, but one can still obtain insights into the physics near the 
CEP. In this contribution we focus on the nature of the soft mode. 



2. The method for solving the model 

In the presence of the vacuum expectation value $ the requirement of a 
finite constituent quark mass m q — g<& as TV — * oo defines the Lagrangian 



L[a,7r a ,V]=- 
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p 2 = a 2 + it 2 and 

/i is the external field. The fcrmions, for which the chemical potential 
is introduced by the shift d t — > 9 t — i/i, start contributing in the large 
N expansion at level 0(1/ y/N), which is between LO and NLO of the 
mesonic sector. To calculate their corrections to the LO, which has O(N) 
symmetry, we disregarded the second independent quartic coupling of the 
model, present for Nf > 2. 

Due to the proliferation of graphs a resummation is needed. The best 
thing to do would be to use a selfconsistent pion and fermion propagator. 
This makes renormalisation difficult, and a consistent method for it was 
proposed only recently 2 . Our method consists in taking into account the 
N c coloured quarks to 0(1/ y/N) perturbatively at one loop order and using 
for the propagator of the pions a simplified form, parameterised by a mass 
M determined from the gap-equation: 

,9 2 N C 



M 2 = -iG-\ P - 



0) = m 2 + ^$ 2 + ^T B (M) 
o 6 



T F (m q ). (2) 



The pion and fermion tadpoles (Tb,f) determine the equation of state (EoS) 



$ff($ 2 )-/i = $ 



m 2 + ^($ 2 + r B (M)) 



4g 2 JV e 



T F (m q ) 



i = 0. 

(3) 

The renormalisation method was presented in 3 . The relation M 2 = H (<f> 2 ) 
shows that in the symmetry broken regime defined by if(<f> 2 ) = and 
in the chiral limit (h = 0) the physical pion mass M is always zero, i.e. 
Goldstone's theorem is fulfilled in the minimum of the effective potential. 
In what follows the number of flavors is set to two: Nf = 2. 
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3. Chiral case (h — 0): analytical determination of TCP 

Since in this case M = 0, one obtains 4 the phase boundary by expanding 
the fermion tadpole in powers of The 2 nd order line is given by 

A <7 2 ^V 

72 



2 2 



4tt 2 



Nr. 
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T = 0, 



(4) 



while the location of TCP comes from combining (4) with the condition 
d 



A , = A g 4 N c 
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n=o Mqb 



= 0. (5) 



lnci/2 = 1 — je + Vi an d ^ = lnMos/MoF gives a relation between M b 
and Mqf the renormalisation scales of the bosonic and fermionic tadpoles. 
The parameters are fixed at T — (i = 0. The choice m q o = toat/3 rs 312.67 
McV and $ = /tt/2 gives g = 6.72. A = 400 is fixed by requiring the 
agreement between the location of the complex cr-pole on the 2 nd Riemann 
sheet and the experimentally favoured mass and width of the a particle. 
Solving (4) and (5) the 2 nd order phase transition, the first spinodal line 
and the location of TCP are found (l.h.s. of Fig. 1). On the r.h.s. we see 




Figure 1. L.h.s.: The T — fi phase diagram for Mob = 886 MeV and r\ = 0. R.h.s. the 
dependence of T c (/i = 0) (upper curves) and of Ttcp (lower curves) on Mqb- 



the variation of TCP and of the critical temperature with rj and M b ■ The 
pole structure of the a propagator restricts the allowed range of M g to 
the right from the arrow, seen on the horizontal axis, where no low scale 
imaginary poles appear besides the well known large scale tachyon, whose 
scale increases with decreasing values of A. Although in this region one can 
easily obtain T c (/z = 0) e (150, 170) McV, the TCP stays robustly below 
70 MeV. Previous effective model studies 5 give similar low values compared 
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with the CEP temperature obtained in 1 , which could mean that the phase 
transition is driven by higher excited hadronic states rather than by the 
light d.o.f. This is supported by a resonance gas model calculation 6 used 
to reproduce lattice results on quark number susceptibilities and pressure. 

An approximate effective potential is obtained by integrating the EoS 
away from equilibrium and subsequently Taylor expanding it around $ = 0. 
This value becomes a minimum as T — > T c and so M — > 0. With this, one 
can show the vanishing of the effective coupling in the large N limit 7 at 
the critical point. Here we just quote the result of 8 



Kpprox($) - 4[/ 2 W 



2 + 12 + \2lQm 2 T + 3 ' 



(6) 



where K = WFWW^ [Li„(-e^) +Li„(-e-*)] |„ = _ 2 , W = 1 - x 
In SjiL and In || = 5—7. The shape of V appmx (<fr) depends on the values of 
the parameters. In the temperature range of a 2 nd order phase transition 
W < 0. For increasing values of fi/T, n changes sign at 1.91. This sign 
change restricts the location of the TCP to the region \i > 2T, where k > 0. 
In the broken symmetry phase < and it vanishes at T c signalling that 

the coefficient of the quartic term in $ also vanishes. Around T c from (4) 

2 

one gets m\ ~ T — T c , which through m e ff = | | gives the correct LO 
critical exponent v = 1 for the O(N) model at large N. The usual Landau 
type analysis applies to the square bracket in (6). The scaling exponent of 
the order parameter on the 2 nd order line is given by the first two terms. 
The minimum condition gives $ 2 = -^i $ ^ (T - T C Y ', = \. At 
TCP one can set A = in (6) and keep the sixth order term. The minimum 
condition gives $ 4 = => $ - (T - T c f , = \. This mean field 

estimates for were checked numerically using the exact EoS (3). 



4. Case of the physical pion mass: the soft mode at CEP 

In view of the gap-equation (2) and EoS (3) the form of the a propagator 

Y±l± 

d* 2 



can be inferred from the consistency condition — G a 1 (p = 0) = d v tl ! S^ \ 



where Ib,f are the pion and fermion one-loop "fish" integrals. Requiring 
M = 140 MeV at T = \i = one has from (2) h/fl = 1.13. The location 
of CEP is determined numerically by simultaneously solving (2) and (3): 
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(/xb,7 1 )cep = (987, 12.34) MeV. Tcep is unrcalistically low, nevertheless it 
is of interest to study which mode becomes soft at this point. Two minima 
meet at CEP which flattens the effective potential. This is shown also by 
the diverging peak of the susceptibility d^/dh. This static information 
has to show up consistently when investigating G a {p) in the static limit 
i.e. first taking p — > 0, then |p| — > 0. The question is in which (p , |p|) 
region one has to look for the responsible excitations. The behaviour of the 
spectral function p(po, |p|) = — ^ lim Im G CT (po + |p|) f° r fixed |p| gives 
a hint in this respect: a peak is found in the po < |p| region. The location 
of its maximum moves toward the origin with decreasing values of |p| and 
at CEP eventually the peak diverges at p n = as |p| — ► 0. The nature 
of the analytical object on the second Riemann sheet which produces this 
peak will be explored further. In the NJL model a pole on the negative 
imaginary axis, approaching the origin as |p| — ► 0, is responsible for the 
peak observed in the spectral function 9 . 

5. Conclusions 

Contrary to the chiral case, where the zero temperature a pole continuously 
moves to the origin of the 2 nd Riemann sheet and becomes soft at TCP, 
for a physical pion mass the a pole stays massive at CEP. A space-like 
(0 < po < |p|) excitation is responsible for the flattening of the effective 
potential at this point. Before a more detailed investigation of the analytic 
structure at CEP we have to improve our method of solving the model 
by dynamically generating the fermion mass through a Dyson-Schwinger 
equation for the fermion propagator, but also one has to extend our method 
to the physically more appealing SU(3)l x SU(3)r case. 
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